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RADAR  PULSE  COMPRESSION  AND 
ELECTROMAGNETIC  INTERFERENCE  (EMI) 


1.  INTRODUCTION 

Electromagnetic  interference  (EMI)  is  a  well-known  problem  that  affects  most  radar  and  tele¬ 
communication  systems.  In  the  past,  radars  have  had  relatively  narrow  bandwidths  and  EMI  between 
radars  was  kept  to  a  minimum  by  frequency  allocation  to  different  frequency  bands.  However,  the 
current  trend  is  toward  wide  operating  agile  and  instantaneous  bandwidths,  possibly  50%  or  more, 
which  necessarily  overlap.  EMI  reduction  therefore  requires  other  means  than  frequency  allocation  to 
different  bands,  possibly  noninterfering  coded  waveforms,  as  pursued  in  this  report. 

EMI,  in  a  broad  sense,  includes  intrasystem  problems  as  well  as  intersystem  problems.  A 
number  of  solutions  have  been  used  or  proposed  to  reduce  EMI,  based  on  the  characteristics  of  radi¬ 
ated  interference,  transmission  path,  and  receiver  configuration.  They  also  depend  on  how  the  prob¬ 
lems  are  defined,  that  is,  from  the  viewpoint  of  signal  transmission,  signal  reception,  or  both.  The 
general  methodology  may  include:  (1)  grounding  and  bonding,  (2)  shielding,  (3)  use  of  radiation 
absorption  material  (RAM)  or  chemical  agent,  (4)  filtering,  (5)  frequency  management,  (6)  interfer¬ 
ence  cancellation  system  (ICS),  (7)  equipment  design  such  as  the  use  of  multicouplers,  (8)  advanced 
signal  processing  via  state-of-the-art  technology  and  expert  systems,  and  many  others.  From  a  radar 
aspect,  the  electromagnetic  interference  is  considered  to  be  externally  generated  and  does  not  include 
clutter  and  intrasystem  noise.  Generally,  differences  between  desired  signals  and  multipath/externally 
generated  signals  can  be  used  to  suppress  the  interference  and  to  augment  the  desired  signals.  Typical 
interference  suppressers  are  sidelobe  blankers,  defruiters,  adaptive  arrays,  mainlobe  notchers, 
transmission/predetection  filters,  coherent  sidelobe  cancellers  or  adaptive  interference  cancellers  [1], 
which  directly  couple  a  sample  of  the  interference  waveform  from  the  transmitter  sources  to  the 
receiver. 

To  date,  besides  the  above  methodology,  not  much  work  has  been  reported  in  eliminating  EMI 
from  the  viewpoint  of  basic  signal  structure.  In  this  report,  we  propose  a  fundamental  approach  to 
eliminating  EMI  by  applying  complementary  coded  signals.  We  consider  transmitted  multiple-pulsed 
(dissimilar)  waveforms  as  a  burst  of  pulses  that  are  pulse-compressed  by  complementary  sequences. 
We  then  search  a  class  of  these  characterized  signals  in  which  the  property  of  mutual  orthogonality  is 
attained.  As  a  result,  different  signals  in  a  specifically  derived  set  can  be  assigned  to  different  radars 
on  the  same  or  different  platforms  in  which  the  signals  are  mutually  orthogonal  in  the  sense  that 
cross-correlation  of  signals  is  zero.  Therefore,  the  interference  from  other  radar  sources  may  be 
reduced.  Since  the  derived/applied  signals  are  all  complementary  coded  waveforms,  the  individual 
signal,  as  we  refer  to  its  composed  multiple  pulses,  is  presented  in  an  orthogonal  matrix  form.  This 
implies  that  the  signal  autocorrelation  (without  considering  Doppler  shift)  is  zero  at  least  over  the 
unambiguous  range  except  at  the  matched  point.  Thus,  zero  sidelobes  of  the  stationary  targets  are 
attained.  Specific  waveforms  in  this  category  have  the  above  property  extended  beyond  the  maximum 
unambiguous  range.  When  the  signal  is  continuously  transmitted  and  matched  to  the  reference  for 
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specific  range  intervals,  then  the  stationary  clutter  or  targets  from  the  mismatched  range  intervals  are 
eliminated.  The  methodology  described  here  provides  a  means  of  fundamentally  eliminating  EMI 
without  degrading  the  radar  performance. 

In  Section  2,  effects  caused  by  the  use  of  pulse  compression  waveforms  against  electromagnetic 
interference  are  described.  We  consider  a  practical  example  of  using  a  transmission  filter  for  elim¬ 
inating  interference  to  the  Air  Traffic  Control  Radar  Beacon  System  (ATCRBS)  airborne  Mode  S 
transponder  and  Traffic  Alert  and  Collision  Avoidance  System  (TCAS)  receiver  in  an  upgraded  SPS- 
49  system.  Then,  in  Section  3,  we  detail  the  fundamental  approach  to  eliminating  EMI  through  the 
use  of  complementary  coded  signals.  Advanced  waveform  synthesis  is  included  in  this  section  to 
solve  imperfection  of  orthogonality  in  the  presence  of  Doppler  shift. 

2.  TRADEOFFS  BETWEEN  RADAR  PERFORMANCE  AND  EMI  SUPPRESSION 

It  is  the  current  trend  in  radar  design  to  have  a  relatively  wide  bandwidth  (both  operating  agile 
bandwidth  and  instantaneous  bandwidth).  As  a  result,  very  high  resolution  can  be  achieved  for  radars 
to  provide  improved  clutter  suppression  for  the  detection  of  small  targets.  Required  peak  power  is 
reduced.  In  addition,  unwanted  EMI  can  be  notched  out  through  filtering  or  ICS.  Pulse  compression 
waveforms  are  thus  widely  used.  The  commonly  used  waveforms  in  radar  are  frequency,  phase,  or 
amplitude  modulated.  The  popular  pulse  compression  waveforms  are: 

(1)  linear  frequency  modulation  (chirp)  waveforms, 

(2)  Barker  and  pseudorandom  binary  phase  codes, 

(3)  step  frequency  modulation  (step  chirp)  waveforms, 

(4)  step-frequency-derived  polyphase  codes  such  as  Frank  and  P\  codes, 

(5)  linear-frequency-derived  polyphase  codes  such  as  and  P4  codes, 

(6)  complementary  codes,  and 

(7)  Huffman  codes. 

In  the  above,  the  complementary  codes  can  be  binary  or  polyphase  codes  and  the  Huffman  codes 
are  amplitude  and  phase  modulated.  For  the  unweighted  linear  chirp  waveform,  its  autocorrelation 
function  has  the  form  of  sin  x  /x,  which  shows  high  peak  sidelobes.  However,  the  sidelobe  level  can 
be  brought  down  by  applying  appropriate  weighting  functions,  e.g.,  -42.8  dB  peak  sidelobe  by  using 
a  Hamming  weighting.  The  chirp  waveforms  have  range-Doppler  coupling  properties  and  are 
Doppler  tolerant. 

The  binary  codes,  however,  have  relatively  poor  sidelobes.  For  an  aperiodic  pseudorandom 
code  of  length  N,  the  peak  voltage  sidelobe  best  obtainable  is  about  1/A.  It  is  the  ideal  case  that  the 
voltage  sidelobe  level  reaches  \/N2  as  seen  in  the  Barker  or  periodic  pseudorandom  codes.  However 
the  binary  phase-modulated  waveforms  are  not  Doppler  tolerant.  Other  codes  listed  in  (3)  through  (7) 
possess  specific  characteristics  and  may  show  combined  favorable  properties  different  from  the 
forementioned  codes.  For  example,  the  polyphase  Frank  codes  have  lower  voltage  sidelobes  of  about 
1  /f 7r  N)  and  are  more  Doppler  tolerant  than  the  binary  codes.  Many  useful  radar  waveforms,  with 
associated  properties,  are  described  in  the  literatures  (e.g.  [2-4 J ) . 
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In  Section  3,  we  discuss  in  detail  the  complementary  codes,  which  when  considered  as  multiple- 
pulsed  (dissimilar)  signals,  appear  to  be  promising  for  both  suppressing  EMI  and  detecting  a  target. 
In  this  section,  the  EMI  solution  is  restricted  to  the  conventional  approach  through  the  use  of  a 
transmission  filter.  Actually  the  use  of  a  desired  filter  to  eliminate  interference  will  produce  other 
effects  detrimental  to  the  radar  performance.  The  portion  of  the  signal  that  is  filtered  away  is  no 
longer  part  of  the  energy  available  for  target  detection.  Also,  a  distortion  of  the  pulse  compression 
waveforms  being  used  is  unavoidable.  These  effects  cause  power  loss  and  internal  subcode  interfer¬ 
ence.  Here  we  examine  how  the  radar  performance,  i.e.,  the  matched  filter  output,  is  affected  by  the 
application  of  an  EMI-suppressing  filter  and  pulse  compression  waveforms. 

In  conjunction  with  the  AN/SPS-49  waveforms  and  system  architecture  study  [5],  we  assume 
that  the  operation  bandwidth  is  850  to  1400  MHz  for  a  possible  upgraded  system.  We  also  assume 
that  a  360  kW  transmitter  is  used.  According  to  Ref.  6,  the  victim  receivers  considered  for  interfer¬ 
ence  are  the  ATCRBS  airborne  Mode  S  transponder  and  the  TCAS  receiver.  To  eliminate  EMI  as 
seen  from  the  aspect  of  signal  transmission,  a  transmission  filter  (instead  of  a  predetection  filter)  with 
the  following  characteristics  must  be  applied.  (See  Table  1.) 


Table  1  —  Filter  Rejection  Characteristics 


Rejection  Magnitude 
dB 

Rejection  Bandwidth 
MHz 

ATCRBS  Mode  S  Transponder  (1030  MHz) 

55 

44 

TCAS  <  Receiver  (1090  MHz) 

50 

22 

To  examine  how  the  matched  filter  output  is  affected  by  the  this  designed  filter,  we  consider  the 
applied  waveform  to  be  a  step-frequency-derived  polyphase  Frank  code  with  a  code  length  of  100. 
This  polyphase  Frank  code  is  the  concatenation  of  unit  pulses  with  their  phases  equal  to 

<t>ij  =  ~U  —  l)(y  —  1);  y  =  1,2...  10  and  /  =  1.2,.,  10.  (1) 

The  reason  to  consider  this  waveform  is  that  the  Frank  code,  when  applied  as  multiple-pulsed  signals, 
is  also  a  complementary  code  that  appears  to  be  promising  in  suppressing  EMI  (see  Section  3).  It  is 
known  that  the  receiver  matched  filter  output,  in  the  absence  of  noise  and  Doppler  shift,  is  the 
response  of  signal  autocorrelation;  that  is,  the  convolution  of  c(f)  with  c*U)  as  defined  in  Fig.  1(a), 
where  c(t)  and  c*(t)  are  the  time  reversal  and  complex  conjugate  of  r(/),  respectively.  The  auto¬ 
correlation  functions  can  also  be  obtained  from  the  signal  spectra  through  discrete  Fourier  transforma¬ 
tion  (DFT)  as  shown  in  Fig.  1(b). 

The  signal  power  spectrum  of  the  100-element  Frank  code  is  shown  in  Fig.  2(a).  In  the  figure, 
the  coded  signal  is  oversampled  five  times  and  padded  with  499  zeroes  and  then  DFT  processed 
according  to  Fig.  1(b).  Therefore,  there  are  999  samples.  The  abscissa  of  the  figure  indicates  the 
index  over  the  frequency,  and  the  frequency  increment  is  the  reciprocal  of  the  sampling  number  (999) 
multiplied  by  the  sampling  spacing  in  time.  The  envelope  of  this  waveform  spectrum  has  the  form  of 
(sin.r/.tr,  and  the  first  null  appears  at  the  frequency  of  l/x  where  t  is  the  designed  compressed- 
pulsewidth.  The  signal  spectrum  is  further  inverse-DFT  processed  to  obtain  the  signal  autocorrelation 
function,  or  the  matched  filter  output,  as  shown  in  Fig.  2(b).  In  the  figure  and  the  corresponding 
ones  through  Fig.  7(b),  the  compressed  pulse  responses  are  obtained  by  sampling  the  wavform  once- 
per-code-element  (or  per  reciprocal  bandwidth).  The  peak  sidelobe  is  about  -30  dB. 
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Fig.  1(b)  —  An  autocorrelation  function  derived 
through  DFT  processing  of  signal  spectra 


As  noted  before,  we  assume  that  the  signal  spectrum  of  Fig.  2(a)  has  a  bandwidth  from  850  to 
1400  MHz,  approximately  between  the  mainlobe  nulls.  Then  the  interference  located  at  1030  MHz 
and  1090  MHz  (Table  1)  can  be  easily  visualized.  We  consider  first  that  the  interference  is  notched 
out  by  narrow  bandstop  filters  of  1  MHz  centered  at  1030  MHz  and  1090  MHz.  Fig.  3(a)  illustrates 
the  resulting  transmitted  signal  spectrum.  Following  Fig.  1(b)  and  performing  DFT  processes,  we 
obtain  a  matched  filter  output  as  shown  in  Fig  3(b).  The  figure  indicates  nearly  the  same  sidelobe 
level  as  in  Fig.  2(b)  and  shows  little  performance  degradation  because  of  the  use  of  a  cancellation 
filter. 


If  the  interference  needs  to  be  notched  out  by  a  filter  of  larger  bandwidth  such  as  specified  in 
Table  1 ,  then  the  corresponding  transmitted  signal  and  the  resulting  matched  filter  output  are  obtained 
in  Fig.  4.  It  is  shown  that  the  peak  sidelobe  of  the  output  response  is  raised  to  approximately 
-19  dB.  Apparently  the  radar  detection  capability  is  substantially  affected  by  the  insertion  of  a  filter 
for  suppressing  EMI.  In  a  practical  application,  we  may  compensate  somewhat  by  using  a  signal  hav¬ 
ing  a  longer  code  length.  Figures  5  through  7  (corresponding  to  Figs.  2  through  4)  show  the  power 
spectra  and  the  associated  matched  filter  outputs  for  the  Frank  coded  signal  with  a  code  length  of  196 
elements.  The  peak  sidelobes  are  -32  and  -  24  dB,  respectively,  in  the  above  narrow-  and  wide- 
passband  notching  cases  (Fig.  6  and  7). 
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200  400  600  800 

INDEX  OVER  FREQUENCY 

Fig.  3(a)  —  100-element  Frank  code  spectrum 
with  very  narrow  frequency  notching 


SAMPLE  NUMBER  (CODE  ELEMENT) 

Fig.  3(b)  —  The  matched  filter  output  for  the 
Frank  code  spectrum  of  Fig.  3(a) 
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Fig.  4(a)  —  100-element  Frank  code  spectrum  with 
very  wide  frequency  notching 


Fig.  4(b)  —  The  matched  filter  output  for  the 
Frank  code  spectrum  of  Fig.  4(a) 
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Fig.  5(a)  —  1%'element  Frank  code  spectrum 


o  ' 

o 


o 


SAMPLE  NUMBER  (CODE  ELEMENT) 


Fig.  5(b)  —  The  matched  filter  output  for  the 
Frank  code  spectrum  of  Fig.  5(a) 
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Fig.  6(a)  —  196-element  Frank  code  spectrum  with 
very  narrow  frequency  notching 


o 

o 


o 


SAMPLE  NUMBER  (CODE  ELEMENT) 

Fig.  6(b)  —  The  matched  Filter  output  for  the 
Frank  code  spectrum  of  Fig.  6(a) 
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Fig.  7(a)  —  1%-element  Frank  code  spectrum  with 
very  wide  frequency  notching 


Fig.  7(b)  —  The  matched  Filter  ouput  for  the 
Frank  code  spectrum  of  Fig.  7(a) 
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3.  A  FUNDAMENTAL  APPROACH  TO  ELIMINATING  EMI 

Pulse  compression  waveforms  are  often  required  in  a  modern  radar  design.  It  may,  however, 
degrade  the  radar  performance  if  a  conventional  methodology  is  used  for  eliminating  EMI  such  as  dis¬ 
cussed  in  Section  2.  Among  the  pulse  compression  waveforms,  the  signals  derived  from  orthogonal 
matrices  show  potential  applications  in  eliminating  the  multiple  time-around  clutter  [4,  7],  As  to  solv¬ 
ing  EMI  problems,  we  further  examine  these  characteristic  waveforms  and  propose  a  fundamental 
approach  to  eliminating  EMI.  We  search  a  large  class  of  signals  derived  from  orthogonal  matrices  in 
which  the  property  of  mutual  orthogonality  is  attained.  In  doing  so,  different  signals  derived  in  a 
specific  set  can  be  practically  assigned  to  different  radars  on  the  same  or  different  platforms.  A  radar 
alone  processes  the  signals  and  suppresses  the  ambiguous  range  targets  or  clutter.  Among  different 
radars,  interference  is  eliminated  at  least  over  the  unambiguous  range  due  to  mutual  orthogonality  of 
the  multiple-pulsed  waveforms  being  used. 

3.1  A  Class  of  Complementary  Codes 

A  measure  of  pulse  compression  waveforms  to  resolve  signals  in  the  spectral  and  temporal 
domains  is  provided  by  the  receiver  matched  filter  output  or  the  signal  ambiguity  function, 

*> 

X(r,  oj)  =  j.t *(t)  x(t +t)  e~ju>'(lt  (2) 

For  convenience,  we  define  in  (2)  that  x(r,w)  =  2.  In  the  absence  of  noise  and  Doppler 

shift  (u  =  0),  the  above  function  \p(r,0)  is  actually  the  signal  autocorrelation  function.  It  was  shown 
in  [8],  for  two  c<  <  iplementary  sequences  such  as  the  following  pairs  of  Golay  sequences  of  order  16, 

€>=(10  0  01  101  10000010] 

(3) 

b  =  [0  100000101001  1  10] 
that  the  sum  of  their  autocorrelation  functions  has  zero  sidelobes. 

Generally  the  complementary  codes  are  defined  to  be  a  complementary  set  of  sequences  ]9], 
which  can  also  be  represented  by  an  array  matrix.  For  example,  referring  to  the  Frank  codes 
described  in  Eq.  (1),  we  can  form  a  multiple-pulsed  waveform  that  consists  of  i  ( i  =  10)  pulses  with 
each  pulse  compressed  by  the  coded  sequences  {exp  j  =  1,2,...,  10]  for  /  =  1.2..,  10.  Conse¬ 

quently,  the  array  matrix  characterizing  this  multiple-pulsed  waveform  has  the  form  of  \a,j\  with 
Ujj  =  exp  (<t>ij)  indicating  the  respective  elements  at  the  /th  row  and  jlh  column.  Incidently,  the 
above  polyphase  waveform  (or  the  array  sequences)  is  complementary  coded  and,  hence,  the  sum  of 
the  autocorrelation  functions  of  the  above  individual  pulses  is  zero  except  at  the  matched  point. 

For  clarity,  we  consider  another  example  with  a  binary  signal  represented  by  the  following  com¬ 
plementary  code  (in  an  array  matrix  form). 
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Figure  8  shows  the  aperiodic  autocorrelation  function  of  the  above  signal.  It  can  be  seen,  to  the  max¬ 
imum  unambiguous  range  interval,  that  the  sum  of  the  autocorrelation  functions  of  all  array  sequences 
is  zero  except  at  the  matched  point.  Note  that  the  correlation  functions  appearing  in  the  second,  third, 
.  .  time-around  intervals  are  generally  nonzeros. 
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Fig.  8  —  A  periodic  autocorrelation  of  the  multiple  signals  represented 
by  the  complementary  array  code  A 


A  matrix  to  be  sufficiently  a  complementary  array  must  be  an  orthogonal  matrix  in  a  sense  that 
the  column  vectors  of  the  matrix  are  mutually  orthogonal  [4],  Consequently,  the  Hadamard  matrices, 
the  Frank  and  P4  code  square  matrices,  and  matrices  generated  from  a  pseudorandom  shift  register 
code  (such  as  A  above)  contain  complementary  sequences.  Matrix  A  is  formed  after  permutation  by  a 
7  x  7  submatrix  consisting  of  a  simple  pseudorandom  code  (1  —I  —1  —1  1  1  —  1 )  and  its  rota¬ 
tions  and  an  eighth  row  and  column  of  Is. 


It  is  interesting,  as  noted  in  Ref.  4,  that  the  polyphase  Frank  code  and  P4  code  square  matrices, 
when  considered  as  multiple-pulsed  waveforms,  have  zero  autocorrelation  (in  a  periodic  sense)  beyond 
the  maximum  unambiguous  range.  Therefore  when  the  signal  is  continuously  transmitted  (in  a 
periodic  sense)  and  matched  to  the  reference  for  a  specific  range  interval,  then  the  stationary  clutter 
from  the  mismatched  range  interval  is  eliminated.  If  the  filter  is  matched  to  the  most  recently 
transmitted  code,  the  sidelobes  of  stationary  targets  in  the  unambiguous  range  are  zeros.  In  a  sense  of 
zero  cross-correlation,  instead  of  zero  autocorrelation,  the  above  complementary  coded  waveforms 
may  be  applied  to  a  multisensor  operation  so  that  EMI  is  eliminated.  However,  a  search  of  “orthog¬ 
onal  sets”  containing  similar  complementary-coded  signals  is  needed.  Consequently,  different  signals 
in  a  derived  set  can  be  assigned  to  different  radars  in  which  the  applied  signals  are  mutually  orthogo¬ 
nal.  Interference  is  therefore  eliminated  among  different  radar  sources. 

An  orthogonal  set  5  may  be  formed  through  orthogonal  transformation  of  a  specific  array  matrix 
of  interest.  For  example,  let  the  following  An  be  a  k  x  n  complementary  array  matrix  (here, 
k  =  n  =  4)  and  5  =  \A  , ,  A2,  .  .  1  for  which  A,,  A  2,  ■  ■  are  the  matrices  to  be  searched  for: 
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1  1  1 
1  -1  1 
1  -1  1 
1  1  1 

Then  one  possible  5  (corresponding  to  a  Walsh  vector  set  W)  can  be  obtained  by  A,  -  A„  7}  with  T, 
an  n  x  n  diagonal  matrix  and  7}/  —  V,  e  W,  where  V,  is  a  Walsh  vector  and  /  is  an  n-dimensional 
identity  vector.  The  procedures,  with  A„  given  above,  are  illustrated  as  follows.  For  a  Walsh  vector 
set  defined  by  W  =  { V , ,  V2,  F3,  V4\  with  V,  =  (1  1  1  l)r,  V2  =  (1  1  -1  -l)r. 
=(1  -1  1  -l)7,andK4  =(  — 1  1  1  -  1  )r,  we  can  easily  derive  that 

-1-1  11 
1-1-11 
-1  1-11 
till 


-1  1  1-1 
1  1-1-1 
-1  -1  -1  -1 
1-1  1-1 

and 

1-1  1-1 
-1  -1  -1  -1 
1  1-1-1 
-1  1  1-1 

In  the  above,  5  =  \A  \ ,  A2.  A),  A4\  and  the  multiple-pulsed  waveforms  represented  by  the 
matrices  A  | ,  A2.  A2  and  A4  are  mutually  orthogonal.  In  other  words,  any  pair  of  the  signals  with 
their  representing  matrices  belonging  to  S  have  zero  cross-correlation  functions.  Moreover,  each  sig¬ 
nal  is  represented  by  an  orthogonal  matrix  so  that  its  autocorrelation  is  zero  except  at  the  matched 
point.  Figure  9  is  the  aperiodic  autocorrelation  function  for  the  signal  represented  by  /1:.  Figure  10 
illustrates  aperiodic  cross-correlation  between  signals  represented  by  A2  and  A , . 

In  the  example  shown  above,  the  derived  set  S  is  noj  unique.  For  instance,  let  the  Walsh  vector 
be  W  =  IF,,  V2,  Vx,  V4\,  or  W  =  \V , ,  V2,  Vy.  VA\  with  F,=|-l  -1  -1  -l|r. 
V2  -  [  - 1  -I  I  1|;,  Vy=\-\  -1  -I  l]r.  and  K4  =  |1  -1  -1  l)r.  Then  two  more  sets 
containing  mutually  orthogonal  signals  are  obtained.  In  fact,  a  total  of  32  sets  can  be  found  in  this 
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example.  The  above  procedures  can  be  applied  to  a  different  seed  matrix  Aa  to  obtain  a  new  class  of 
multiple-pulsed  signals  of  interest.  If  A(,  is  a  polyphase  Frank  code  matrix,  additional  properties, 
such  as  the  periodic  autocorrelation  functions  being  zero  over  ambiguous  range  intervals,  are 
observed. 


3.2  Advanced  Waveform  Synthesis 

As  discussed  in  previous  sections,  mutually  orthogonal  waveforms  are  fundamentally  applicable 
in  solving  EMI  problems.  With  appropriate  use  of  multiple-pulsed  signals  selected  from  a  derived 
orthogonal  set  containing  complementary  codes,  the  radar  could  suppress  interference  from  other 
radar  sources  and  possibly  eliminate  ambiguous  range  targets  or  clutter.  However,  orthogonality  of 
the  complementary  codes  are  not  preserved  in  the  presence  of  Doppler  shift.  That  is,  the  temporal 
sidelobes  of  matched  filter  outputs  are  no  longer  zeros.  The  compressed  pulse  peak  is  also  reduced. 
For  example,  if  the  Doppler  shift  in  angle  is  7r  for  the  separation  between  two  complementary 
sequences,  then  the  sidelobes  of  the  autocorrelation  of  these  two  sequences  will  add  up  instead  of  sub¬ 
tract  from  each  other.  To  minimize  this  problem,  we  examine  the  option  of  using  orthogonal 
multiple-pulsed  waveforms  incorporated  with  an  inner-outer  coding  scheme.  Here  a  train  of  pulses, 
called  the  outer  code,  is  “modulating”  previously  discussed  complementary  codes  (the  inner  code). 
The  inner  code  sequences  are  indeed  “determined”  by  a  preselected  outer  code  as  described  next. 
The  outer  code  may  also  be  called  a  macro-sequence  [10]. 

A  synthesized  pulse  train  is  simply  shown  in  Fig.  1 1 .  The  pulse  train  consists  of  N2  pulses  with 
a  pulse  repetition  interval  T.  Each  pulse  is  pulse-compressed  by  a  complementary  sequence  confined 
to  a  generalized  array  matrix.  The  selection  of  which  complementary  sequence  to  be  modulated  is 
determined  by  a  predefined  outer  code.  Assume  each  complementary  sequence  has  A,  subpulses  with 
a  subpulse  length  of  5  (see  Fig.  1 1).  Then,  for  the  &th  pulse. 


N  |-1 

«*  =  E  MO*  . 

n  =0 

where  h„(t)  =  1  and  6n  =  0  or  ir  with  n5  <  t  <  (n  +  1)5  and  n  =  0,  1,  .  .  ,  N\  —  1.  We  denote 

n2  - 1 

the  resulting  waveform  by  £  uk(t)  with  kT  <  t  <  kT  +  N  j6  and  k  =0,  1,  .  ,  N2  —  1. 

k  =0 

According  to  Eq.  (2),  the  matched  filter  outputs  are  related  to 


or 


n2t 

\p(T,  W)  =  j 
0 


N2~  1 

E  «*(/)  uk(t  +  T)eJU'dt 

k=  0 


(9) 


n2- i  KT+N'& 

\ P(r,  o>)  =  E  I  -  kT)  ukU  -  kT  +  T)eJU,dr. 

k  =0  kT 

By  replacing  the  parameter  t  -  kT  by  t,  the  above  equation  becomes 

AN- 1  N*s 

Mr,  a>)  =  E  |  uk(t)uk(t  +  T)ej“'e]“kT  dt 
k  =  0  0 


(10) 
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n2t- 


u0(,> 


U1  (t) 


N2-1 


Fig.  1 1  —  A  synthesized  pulse  train 


or 

N2-  1 

\p(r,  to)  =  £  Uke>UT  (11) 

k  =0 

where 


A',5 

Uk  =  j  U*k(t)  Uk(t  +  T)eju>‘dt.  (12) 

o 

The  results  of  Eqs.  (11)  and  (12)  are  particularly  interesting  since  their  explicit  forms  can  be 
easily  analyzed.  As  noted  above,  the  assignment  of  uk  is  determined  by  a  predefined  outer  code  mk. 
To  be  general,  we  consider  the  complementary  array  matrix  has  the  order  of  4  x  N{  and  is 
represented  by  an  array  matrix  (a(t)  b(t)  c(t)  d(t))T.  Consequently,  one  of  the  four  complemen- 
tary  sequences  would  be  assigned  to  uk  where  k  =  1,2,  .  .  ,  N2-  Naturally,  a  good  choice  of  the 
outer  code  mk  is  the  Quadraphase  sequence  that  can  be  generated  by  the  Barker  code  or  any  other 
binary-phase  codes  [11]  if  a  longer  code  length  is  desirable.  Here  uk  is  assigned  as  follows: 
uk  =  aft).  b(t),  c(t)  or  dft)  if  mk  =  1,-1,  j  or  -j,  respectively. 

As  a  result,  Eq.  (12)  can  be  rewritten  as 

r  (m2k  +  1  )(mk  +  1)  (ml  +  1  )(mk  -  1)  ,  x 

- - - a  ( t )  +  - - - b  ( t ) 

4  -4 


(ml  -  l)(m^  +  j)  .  _  (m2k  -  l)(mk  -  j) 

— c  ( t )  +  - — - d  (0 


-4 j 


4/ 


(m2k  +  l)(mk  +  l)  (m2k  +  \)(mk  -  1) 

- - - a(t+r)  +  - - - b(t  +  t) 

4  -4 


(ml  -  1  )(mk  +  j)  (ml  -  1  )(mk  -  j)  ^  n 

+  - - - C(,  +  t)  +  - — - d(t  +  t) 


-4 j 


4 j 


eju  'dt 


(13) 
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If  we  consider  the  matched  filter  outputs  in  the  unambiguous  range  interval  only,  then 


co)  =  £ 

k  =o 


N-~'  <  (ml  +  1  )2(mk  +  l)2  .  ,  (ml  +  l)2(m*  -  l)2 


16 


-\j/a(T,  CO)  + 


16 


tbir,  co) 


(ml  -  l)2(mA.  +  j)2  ,  ,  x  .  (ml  +  i)2(mk-j)2  (  / 

+  - — - ^c(t,  co)  +  - 77 - —Vj(t,  co) 


16 


-16 


eJaTi,  (14) 


where  ^(r.co)  =  ja*(t)a(t  +  T)eJ“'dt.  The  same  definitions  apply  to  \ph(T,w)y  ^,(7,co)  and 

ij(r,u). 

After  simplification,  one  can  easily  derive 


1 


N-,  -  I 


Mr,  co)  =  —  (l pa(r,  co)  +  \ph(T ,  co)  +  I pt.(r,  to)  +  ^(t,  to))  £  e 

4  it  =0 


jujkT 


1  ^2-1 

+  7  (lM*.  co)  +  ^(r,  co)  -  t/v(7,  co)  -  ^(t,  to))  £  ml  ejUT 
4  k  =o 


A/2 -l 


N,  -  I 


+  -L  (^,(7,  co)  -  ^(7,  co))  £  w*  -  -  (tA,(7,  co)  -  ^(7,  co))  £  jmke 


juikT 


k  =0 


it  =0 


N2  -  I  .  A/2-1  . 

+  —  OMt,  u)  -  1M7,  go))  *£  mkejUT  +  —  (&-(t,  to)  -  ^(7,  co))  £  jmkeJukT  (15) 

4  R  =0  4 


*  =0 


and 


1 

to)  <  — 

A/ 2  -  1 

(\l/„( 7,  CO)  +  \ph( 7,  CO)  +  lA<(T,  CO)  +  lA</(7,  CO)) 

£  eiukT 

1  4 

k  =0 

+  i  J 

A/-.  -  l 

(1/^(7,  co)  +  ypb(7*  co)  -  ^ ,  (t ,  co)  -  ^(7,  co)) 

4  1 

*  =0 

V,  -  1 

(\pa(T,  U))-\j/h(T.  C0)-l/'((7,  C0)  +  l/'(/(7,  CO)) 

*  =0 

A/.  -  1 

j 

(\J/tl(T,  co)  -  'Ph(T'  4  )MT'  w)  _  'PA7-:  GO)) 

*  =0 

(16) 


In  Eq.  (16),  the  summation  terms  are  in  DFT  format,  and  hence  \(/(t, co)  can  be  easily  calculated  if  the 
complementary  array  matrices  and  outer  codes  \mk\  are  given. 
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For  this  specific  case,  with  the  analytical  results  shown  in  Eqs.  (15)  through  (20),  the  matched  filter 
response  can  easily  be  realized. 

To  further  evaluate  these  synthesized  waveforms  and  compare  the  results  consistently,  we  deter¬ 
mine  the  matched  filter  outputs  at  the  presence  of  Doppler  shift  for  the  following  cases  that  the  outer 
codes  are:  (a)  simply  two  pulses,  (b)  a  special  case  of  the  Quadraphase  sequence  generated  from  a 
13-element  Barker  code,  (c)  the  derived  sequence  of  case  (b)  plus  —j,  and  (d)  a  complementary 
sequence  (generated  by  a  pseudorandom  sequence  plus  one).  In  all  cases,  we  select  the  complemen¬ 
tary  array  matrix  An  to  be  the  Golay  mates  of  order  16,  that  is  A0  =  (a  b)T  with  a  and  b  defined  in 
Eq.  (3).  These  two  complementary  sequences  (the  inner  code)  are  then  modulated  by  the  outer  code, 
which  has  a  pulse  repetition  period  of  T  and  a  duty  cycle  of  50%.  A  computer  program  closely  fol¬ 
lowing  Eq.  (11)  is  used  to  generate  the  matched  filter  outputs  for  the  above  cases.  The  Doppler 
phase  shifts  <t>  are  referenced  to  the  entire  code  length  N2T,  i.e.,  <t>  =  2-irf(N2T)  and  is  varied  from  0 
to  approximately  7r. 

Case  A  :  This  is  the  baseline  case  where  the  pulse  train  contains  a(t)  and  b(t)  only.  In  this 
case,  the  matched  filter  outputs  are  very  much  affected  by  the  Doppler  shift.  According  to  Eq.  (1 1), 

\p(T,  z)  =  U0(t,  Z)  +  U\(t,  z)  ej2*z.  (21) 

Clearly.  \p{r,z)  =  U0(t,z)  +  V\ (r,z)  or  ^(r,z)  =  U0(t,z)  -  U i(r,z)  when  z  =  0  or  1/2.  Figure 
12  shows  the  3-D  matched-filter  output  response  \iJ/(t,z)  1 .  In  the  figures,  there  are  31  (or 
2  x  16-1)  samples  in  the  time  axis  and  51  samples  in  the  frequency  axis.  The  magnitude  is  nor¬ 
malized  to  the  response  peak  at  r  =  z  =0,  i.e.,  xmax  =  N^Nl- 
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i 

Fig.  12  —  Ambiguity  diagram  of  the  Golay  complementary 
sequences  of  length  16 


Case  B  :  Here  we  consider  the  outer  code  to  be  a  special  case  of  the  Quadraphase  sequence  gen¬ 
erated  from  a  13-element  Barker  code.  We  assume  the  assigned  complementary  sequences  are 
a(t ),  b(t),  a(t),  and  b(t)  when  the  outer  code  mk  =  1,  -1,  j  and  —j,  respectively.  This  results  in 
i pa(r,  z)  =  z)  and  \pb  (r,  z)  =  \j/d  (r,  z)  in  Eq.  (14),  and  Eq.  (15)  becomes 


Mr,  Z)  <  +  \fb) 

N-,  -  1 

+  ~  tb) 

JVi  -  1 

E  mleMT 

l 

k  =0 

z 

Xr 

II 

O 

With  Eqs.  (19)  and  (20),  Eq.  (22)  can  be  further  simplified  to 


. .  ,  1  , ,  ...  sin  1 37rz 

z)  <  -  (4  +  ib)— - 

2  sin  7rz 


1  , ,  ,  ,  ,  cos  26tz 

+  -zii'a  -  n)  12  +  - 

2  COS  7TZ 


It  is  interesting  to  see  that  the  matched  filter  outputs  are  composed  of  the  sum  and  difference  of  the 
autocorrelation  functions  of  complementary  sequences  weighted  by  wj  =  sin  13irz/sin  7rz  and 
w 2  =  (12  +  cos  267rz /cos  7rz)l/2,  respectively. 


When  t  ^  0,  +  \ph  is  small  due  to  the  property  of  complementary  sequences.  The  filter 

output  result  hence  depends  on  the  second  term  in  which  some  sidelobes  of  \J/a  -  \pb  are  usually  not 
identically  zeros.  Clearly  the  second  term,  and  hence  the  matched  filter  output,  is  minimized  if  we 
can  find  a  good  outer  code  such  that  the  weight  of  \j/a  -  \ph  is  relatively  flat  and  bounded.  It  is 
indeed  the  case  here.  Figure  13  shows  the  ambiguity  diagram  of  the  Golay  sequences  of  length  16 
outercoded  by  a  modified  Quadraphase  code.  It  appears  that  the  output  responses,  except  at  zero 
Doppler  shift,  demonstrate  much  better  performance  compared  with  that  shown  in  Case  A.  For  exam¬ 
ple,  for  a  Doppler  shift  of  0.8ir  over  the  entire  pulse  train,  the  main  peak  sidelobe  is  about  —47  dB 
compared  to  - 19  dB  in  Fig.  12  of  Case  A.  In  both  cases,  the  reduction  of  the  matched  main  peak  is 
nearly  the  same. 
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Fig.  13  —  Ambiguity  diagram  of  the  Golay  sequences  outercoded 
by  a  modified  Quadraphase  code 


Note  that  the  weighting  function  vt>|  is  oscillatory  having  nulls  at  the  frequencies  of  z  =  n/13, 
n  =  1.  2,  .  ..  Even  at  the  matched  point  t  =  0,  the  result  of  t  =  0  is  expected  at  the  above  fre¬ 
quencies.  Also,  in  the  above,  the  resulting  synthesized  waveform  is  actually  a  binary  sequence  (due  to 
simplification)  and  \mk\  =  {a,  a,  b ,  b ,  a,  b ,  a ,  b ,  a,  b,  b,  a,  a\.  There  are  7  a’s  and  6  b’s  in  the 
sequence.  The  unequal  numbers  of  a  and  b  result  in  nonzero  sidelobes  of  the  matched  filter  output  at 
zero  Doppler  shift  (Fig.  13). 

Case  C  :  In  this  case,  we  form  a  new  outer  code,  i.e.,  1,  j,  -1,  -j,  1,  -j,  1,  -j,  1,  -j,  -1, 
j<  ~j<  which  is  obtained  by  adding  —j  to  the  previous  Quadraphase  sequence  generated  from  a  13- 
element  Barker  code.  In  other  words,  one  more  pulse  coded  by  b  is  concatenated  to  the  binary 
sequence  |m*)  derived  in  Case  B.  The  resulting  sequence  contains  equal  numbers  of  o’s  and  b's. 
Zero  sidelobe  response  is  thus  preserved  for  the  matched  filter  output  when  the  Doppler  shift  is 
absent.  Figure  14  shows  the  ambiguity  diagram  of  this  synthesized  waveform.  Clearly  the  matched 
filter  outputs  for  this  synthesized  waveform  perform  superior  to  the  corresponding  ones  as  seen  in 
Cases  A  and  B. 

Case  D  :  The  case  that  the  outer  code  be  a  complementary  as  suggested  in  Ref.  1 1  is  simulated 
here.  We  select  the  outer  code  5  as  a  pseudorandom  sequence  plus  one.  For  the  code  length  to  be 
closer  to  14  as  used  in  the  previous  examples,  we  consider  that  the  pseudorandom  sequence  has  a 
code  length  of  15  and  is  generated  by  a  polynomial  of  023  in  octal.  We  also  choose  the  initial  condi¬ 
tion  of  7  in  octal,  which  results  in  a  lowest  peak  sidelobe  of  3  or  of  10  log  (3/1 5)2  dB  as  normalized 
to  the  mainlobe  peak.  It  is  easy  to  derive  that  S  =  {0  I  I  1  I  0  1  0  1  10  0  10  0  0). 
With  this  outer  code  modulating  the  complementary  sequences  a  and  b ,  the  matched  filter  output 
response  is  shown  in  Fig.  15.  As  before,  we  vary  the  Doppler  shift  over  the  entire  pulse  train  from  0 
to  7r.  The  results  show  improvement  in  radar  performance  but  are  not  as  good  as  in  Case  C  when 
both  are  compared  to  the  baseline  Case  A. 
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Fig.  15  —  Ambiguity  diagram  of  the  Golay  sequences  outercoded 
by  a  16-element  complementary  sequence 
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4.  DISCUSSION  AND  CONCLUSIONS 

In  this  report,  use  of  radar  pulse  compression  against  EMI  is  analyzed.  We  show  that  mutually 
orthogonal  waveforms  specifically  derived  from  complementary  coded  signals  are  fundamentally 
applicable  in  solving  EMI  problems  without  degrading  the  radar  performance.  We  consider  transmit¬ 
ted  multiple-pulsed  (dissimilar)  waveforms  as  a  burst  of  pulses  that  are  pulse-compressed  by  comple¬ 
mentary  sequences.  We  then  search  a  class  of  these  characterized  signals  in  which  the  property  of 
mutual  orthogonality  is  attained.  As  a  result,  different  signals  in  a  specifically  derived  set  can  be 
assigned  to  different  radars  on  the  same  or  different  platforms  in  which  the  signals  are  mutually 
orthogonal  in  the  sense  that  cross-correlation  of  signals  is  zero.  A  radar  alone  processes  the  signals 
and  suppresses  the  ambiguous  range  targets  or  clutter.  Among  different  radars,  interference  is  elim¬ 
inated  at  least  over  the  unambiguous  range  due  to  mutual  orthogonality  of  the  multiple-pulsed 
waveforms  being  used.  Synchronization  among  multiple  radars  may  be  effective  in  reducing  prob¬ 
lems  associated  with  cross-correlation  in  different  unambiguous  range  intervals. 

However,  imperfection  of  orthogonality  is  observed  in  the  presence  of  Doppler  shift.  In  consid¬ 
ering  severe  Doppler  shift,  Doppler  filter  bank  receivers  can  be  practically  implemented  to  process 
the  target  signals.  Otherwise,  this  undesirable  effect  may  be  partially  compensated  for  through 
advanced  waveform  synthesis.  We  describe  the  option  of  using  complementary  coded  signals  incor¬ 
porated  with  an  inner-outer  coding  scheme.  We  generate  the  characteristic  ambiguity  diagram  of  a 
synthesized  waveform  composed  of  a  simple  pair  of  Golay  sequences  outercoded  by  a  pseudorandom- 
code-generated  complementary  sequence.  If  the  outer  code  is  chosen  from  a  modified  Quadraphase 
code,  even  better  waveform  characteristics  can  be  achieved.  Both  synthesized  waveforms  perform 
superior  to  the  baseline  case  where  the  signal  is  simply  two  pulses  and  coded  by  the  Golay  sequences. 

In  conclusion,  with  appropriate  use  of  the  multiple-pulsed  waveforms  selected  from  a  derived  set 
of  complementary  coded  signals,  the  radar  would  suppress  interference  from  other  radar  sources.  In 
addition,  stationary  targets  or  clutter  from  ambiguous/mismatched  range  intervals  may  also  be  elim¬ 
inated.  In  the  future,  we  plan  to  generalize  this  fundamental  principle,  perhaps  through  waveform 
synthesis  and  also  to  search  for  best  seed  matrices  and  various  orthogonal  transformations  in  obtaining 
new  signal  sets  of  interest.  It  would  be  ideal  if  the  described  zero  autocorrelation  (except  at  the 
matched  point)  and  zero  cross-correlation  between  signals  could  be  extended  over  ambiguous  range 
intervals  and  also  be  valid  in  the  presence  of  Doppler  shift. 
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